Steady-State Entanglement for Distant Atoms by Dissipation in Coupled Cavities 
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We propose a scheme for the generation of entangled states for two atoms trapped in separate cav- 
ities coupled to each other. The scheme is based on the competition between the unitary dynamics 
induced by the classical fields and the collective decays induced by the dissipation of two delocalized 
field modes. Under certain conditions, the symmetric or asymmetric entangled state is produced 
in the steady state. The analytical result shows that the distributed steady entanglement can be 
achieved with high fidelity independent of the initial state, and is robust against parameter fluctua- 
tions. We also find out that the linear scaling of entanglement fidelity has a quadratic improvement 
compared to distributed entangled state preparation protocols based on unitary dynamics. 

PACS numbers: 03.67.Bg, 42.50.Pq, 03.67.-a 



There have been various practical applications for 
quantum entangled states, ranging from quantum tele- 
portation P, Q to universal quantum computation [3, 0] ■ 
The main obstacle in preserving entanglement is decoher- 
ence induced by the environment. Recently, dissipative 
state preparation has become a focus in quan tum compu- 
tation and entanglement engineering |5|-|2(]||. which uses 
decoherence as a powerful resource without destroying 
the quantum entanglement. These schemes are robust 
against parameter fluctuations, obtain high fidelity en- 
tanglement with arbitrarily initial states, and do not need 
accurate control of the evolution time. Particularly, Kas- 
toryano and Reiter et al. [1, Q proposed a novel scheme 
for dissipative preparation of entanglement for two atoms 
in an optical cavity which gets a qualitative improvement 
in the scaling of the fidelity with optimal cavity parame- 
ters as compared to any state preparation protocol with 
coherent unitary dynamics. However, most of the previ- 
ous theoretical schemes and experiments |2l| concentrate 
on the case in which two atoms are trapped in a single 
cavity. 

For distributed quantum information processing, it is 
a basic requirement to perform state transfer and quan- 
tum gate operation between separate nodes of a quan- 
tum network. To overcome the difficulty of individual 
addressability existing in a single cavity, efforts have 
been devoted to the coupled-cavity models both theo- 
retically [23 - |28| and experimentally [1^. Most works on 
the coupled-cavity system focused on the traditional co- 
herent unitary dynamics, requiring precise timing and 
special initial states. Clark et al. [30J proposed a scheme 
to entangle the internal states of atoms in separate op- 
tical cavities using technique of quantum reservoir engi- 
neering, however the scheme requires a complex atomic 
level configuration. Furthermore, the evolution towards 
the steady state slows down as the entanglement of the 
desired state increases. 

In this paper, we generalize the idea of Refs. Q 



and propose a scheme for producing distributed entan- 
glement for two atoms trapped in coupled cavities. Due 
to the coherent photon hopping between the two cavities, 
the system is mathematically equivalent to that involv- 
ing two atoms collectively interacting with two common 
nondegenerate field modes symmetrically and asymmet- 
rically, respectively. Each delocalized field mode induces 
a collective atomic decay channel. The present scheme 
uses the competition between the transitions induced by 
the microwave fields and the two collective atomic decay 
channels to drive atoms to a symmetric or asymmetric 
entangled state. Analytical and numerical results show 
that the distributed steady entanglement can be obtained 
with high fidelity. The scheme is independent of the ini- 
tial state and robust against parameter fiuctuations. No 
photon detection, or unitary feedback control is required. 
The linear scaling of -F is a quadratic improvement on the 
cooperativity parameter compared to any known 

entangled state preparation protocol for coupled-cavity 
systems [134111: whose optimal value is 1 — F a C~^/^. 

The experimental setup, as shown in Fig. 1, consists of 
two identical A-type atoms each having two ground states 
|0) and and an excited state |2) and trapped in one 
detuned cavity. An off-resonance optical laser with de- 
tuning A drives the transition |0) |2) and a microwave 
field resonantly drives the transition |0) ^ The cav- 
ity mode is coupled to the |1) o |2) transition with the 
detuning A — (5, where 5 is the cavity detuning from two 
photon resonance. We here assume a Om phase difference 
between the microwave fields applied to the two atoms. 
Under the rotating- wave approximation, the Hamiltonian 
of the whole system in the interaction picture reads Hi 



Ha + H„ 



V+ 



where 



Ho = <5(4ai + 4a2) + A(|2)i(2| -I- |2)2(2|) 
+ [.g|2)i(l|ai + (7|2)2(l|a2 + i/.c.] 
+J{a\a2 + aial), 

ff. = ^(e''"|l>i(0| + |l)2(0|)+i/.c.. 
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FIG. 1. (Color online) Experimental setup for dissipative 
preparation of entangled steady-state between two A-type 
atoms trapped in two coupled cavities. The atom in each 
detuned cavity has two ground states |1) and |0), and one 
excited state |2), which is driven by the same off-resonance 
optical laser. The microwave fields applied to the two atoms 
differ by a relative phase of ^m- 



V- = {V-^y , ai is the cavity field operator in cavity i 
(i — 1,2), J is the photon- hopping strength which de- 
scribes cavity and cavity coupling, g is the atom-cavity 
coupling constant, n and I^m represent the classical laser 
driving strength and the microwave driving strength, re- 
spectively. Om — TT (or 0) guarantees a high fidelity for 
asymmetric steady-state \S) — (|01) — |10))/\/2 ( or sym- 
metric steady-state |r) = (|01) -I- |10))/%/^ ). Let us in- 
troduce two delocalized bosonic modes ci and C2, and 
define asymmetric mode ci = (ai — a2)/V^ and symmet- 
ric mode C2 = (ai + 02) /V^, which are linearly related 
to the field modes of two cavities. In terms of the new 
operators, the Hamiltonian Hq can be rewritten as 

Ho - ^[|2)i(l|(ci + C2) + |2)2(1|(C2 - ci) -I- H.c] 

+{S - J)ctci + {5 + J)c+C2 + A ^ |2),(2|.(4) 



The Hamiltonian describes the asymmetric cou- 
pling for the two atoms to the delocalized field mode ci 
and the symmetric coupling to C2. Due to the photon 
hopping these two delocalized field modes are nondegen- 
erate and each induces a collective atomic decay channel. 
The photon decay rate of cavity i (i = 1,2) is denoted 
as Ki and the spontaneous emission rate of the atoms is 
denoted as 7^ [j = 1,2,3,4). Under the condition ki — 
K2 =n, the Lindblad operators associated with the cavity 
decay and atomic spontaneous emission can be expressed 
asL«i =^ci,L-^ =^C2,L^' = VtT |0)i(2|, = 
10)2(21, L^^ = V7i |l>i(2|, = |1)2(2|. We 
assume 71 =72 =73 =74 = 7/2 for simplicity. 

Under the condition of weak classical laser field, we 
can adiabatically eliminate the excited cavity field modes 
and excited states of the atoms when the excited states 
are not initially populated. To tailor the effective decay 
processes to achieve a desired steady-state, we introduce 
an effective operator formalism based on second-order 
perturbation theory js^]. Then the dynamics of 

our coupled cavity system is governed by the effective 




FIG. 2. (Color online) Two effective models for coherent 
and dissipative interactions among states |00), \S) {\T)) and 
1 11), where two microwave fields cause rapid transitions. The 
atoms decay through the cavity from |00) to |r) and \S) with 
effective decay rates Kci 1 and 1, and from \T) and \S) to 
1 11) with the effective decay rates Kci 2 s-nd Kc^ 2; where ^ 
3> 2 and Kc2,i 2> «;c2,2- 7s,i2, 7S,34, 7t,12 and 7t,34 are the 
effective spontaneous emission rates. The loop-like element 
Ox{^^) {X = 00, S, T) represents the square of the coeffi- 
cient in the corresponding term |X)(X| within -ffe// without 
considering Hg. (a) 6m ~ 0. (b) 6m — n. 



Hamiltonian H^ff and effective Lindblad operator L^ff 
1. 



eff 



V^H^], + iH],],)^V+ + Hg, (5) 



(6) 



where Hj^\j is the inverse of the non-Hermitian Hamilto- 
nian Hmh = Hq — ^ "^^{L^Y . The resulting effective 
master equation in Lindblad form is 



+P(^e//)tie//]}, 



(7) 



H^ff = -Re[^R,]\S){S\ ~ Re[^R2]\T){T\ 

~Re[^Ri]mm + Hg, (8) 



(5+ J)2o2 



= — .„ |r)(oo| 



(9) 

\n){T\, 



(10) 



where Re[ ] denotes the real part of the argument, 
9eff = -^,S =S--K,A =A--7, 



Ri = 
R2 - 

i?3 = 



-{6 -J)iS + J) 
6'g^~A'{S'-J){6'+J)' 
-gJ-S'g^ + A'{S' -J){S' +J) 
[g^-A'{5'~J)][g^-A'{6'+J)r 

gJ^S'g^ + A'jS' -J){S' + J) 
[g^^A'{5'-J)][g^-A'{6'+J)Y 
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Jg2 
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{5 + J),D, 



n -i{5 + J) 
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2A 



(11) 



As shown in Fig. 2 (a) and (b), the loop-hke elements 
Ooo(f^^), Ori^"^) and Os(n^) represent the effective- 
Hamiltonian evolution in three triplet states |00), |T) and 
\S) without microwave fields, respectively. For weak op- 
tical driving f2, H^ff — Hg. There exist two effective 
decay channels characterized by L^j^ and L^^j through 
the two delocalized bosonic modes Ci and C2 as compared 
with the case of Ref. Q in which only one decay chan- 
nel is mediated. It is the photon hopping that lifts the 
degeneracy of the two delocalized field modes and leads 
to the two independent decay channels, ^e// indicates 
the effective decay from |00) to \S) at a rate Kci,i and 
from \S) to 1 11) at a rate 2 caused by asymmetric 
ci mode, and L^^j denotes the effective decay from |00) 
to |r) at a rate ^ and from \T) to |11) at a rate ^ 
caused by symmetric C2 mode simultaneously. The decay 
rates Kci 1 ('^£1,2) ^-nd i^c2 1 ('*C2 2) equal to the square of 
the first (second) coefficient in the right side of Eq. (9) 
and Eq. (10), respectively. Set A^,^ — ^^2 = 0; de- 
cays from 15*) to |11) and from \T) to |11) can be both 
largely suppressed. On the other hand, the microwave 
fields drive the transition between the three states |00), 
|r) (IS")) and |11) for Om = 0(7r). The dynamics of the 
full master equation in Fig. 3 (a) and (b) illustrates that 
we can obtain state \S) or \T) of high fidelity, and the 
time needed for reaching the entangled steady-state |T) 
is about two times as large as that of jS*). This is because 
that the optimal ratio of i/'*ci 2 about 2 times as 
large as Kc2 i/'^c2,2- The errors imposed by all possible 
atomic spontaneous emissions should also be taken into 
account. We apply Eq. (6) again to derive four analytic 
expressions of effective spontaneous emissions with the 
other Lindblad operators U*^ , L'^^ , L'*^ and L'^'^ 



Til 

^eff 



7-73 
^eff 



^[^|i?i||00)(00| + ^|i?2|(|r)(T| + |5)(T|) 



'^\R3\{\T){S\ 



iiJL 

2 ^2^2 



\S){S\)], 
i?i|(|T)(00|-l-|^)(00|) 



2V2 



(|i?2||ii)(r| + |i?3||ii)(5|)]. 



(12) 



(13) 



where | • | denotes modulus of the symbol in it, L 



L]}^ and Ljjf - 



taneous emission for \S) state are 



Lgj^. The operators of effective spon- 



= ii}/,5 = V7i;^ir)(^i, 



Til 

^eff,S 
7-73 _ 7-74 
^eff,S - ^eff,S 



(14) 
(15) 
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FIG. 3. (Color online) The populations of four states IS"), 
|T), |00) and |11) versus the dimensionless parameter gt for 
a random initial state. Both curves are plotted for C = 200, 
K — 7/2, Qm = 20/5, n = 5/20 with A, S and J being the 
optimal values for two entangled steady-states, (a) 6m = 0. 
(b) 9m = TT. (c) The fidelity F^s) for steady-state IS) versus 
C, and the coefficient of the finear scaling in Fjg) as a function 
of C with different ratios K/'y is plotted in the inset. 



and the operators of that for |T) state are 

L:}f^^^L:}f^^^^j^r-^\S){Tl (16) 

^ikT = ^lkT = V7?:^iii)(ri, (17) 



where 

7e// 



(^){(gJ)^ + MAS - ^) + 7^^]^] 

(^2- A(5)2(54 + k2^2) 



(18) 



and Js,i=l,2 = 7T,i=l,2 = 7e///16, 7S,i=3,4 = lT,i=3A = 

Jeff/ 8- Then we use the rate equation to evaluate the 
fidelity for the state (j = S oi T) 
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FIG. 4. (Color online) F\s) in the effective two-qubit system 
versus fluctuations of various parameters, (a) F^s) vs if and 
f; (b) F|s> vs and f . 



Pj = KaPoa ~ (Kb + ^ lj,i)Pj^ 



(19) 



where Pj is the probability to be in the state j. The 
first term on the right side of Eq. (19) represents the 
population decaying into the state j with the rate Ka, 
while the other terms express the po|)ulation leaking out 
of the state j with the rate k;, + Suppose Pj 

~ 1 and the probability in each of the other three states 
is nearly PqOi then 



1-F| 



\s) 



(3 



+ 97e//j/[^^2 



Di 
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Bl 

ki 



(20) 



where i^s) l('5'|pss|S') | is the fidelity of state 15). Set- 
ting Jg^ = A(52 - J^) and K(A5-.gV2) ~ 7((52- j2)/2, 
the optimal fidelity of the entanglement can be obtained. 



The effective two-qubit system in the inset of Fig. 3 (c) 
shows that the fidelity scaling of state 15*) is indepen- 
dent of different ratios k/j, then we find out the actual 
constants for maximizing the fidelity as follows 



F\s) w 12. 8C" 



(21) 



The influences of different parameter fluctuations on 
the fidelity F|5^ of entangled state are considered. As 
shown in Fig. 4 (a) and (b), F^s) keeps above 90% even 
5% fiuctuations in these parameters. The preparation 
process of state |T) is similar to that of \S). 

Photonic band gap cavities coupled to atoms or quan- 
tum dots are suitable candidates for realizing the pro- 
posal. Cooperativity of value C ~ 100 has been realized 
[33| . The cavity modes can be coupled via the overlap of 
their evanescent fields or via an optical fiber, and photon 
hopping between two cavities has been observed |34l |. 

In conclusion, we have proposed a scheme for dissi- 
pative preparation of entanglement between two atoms 
that are distributed in two coupled cavities. We find the 
linear scaling of the fidelity is a quadratic improvement 
compared with distributed entangled state preparation 
protocols based on unitary dynamics. 
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